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ABSTRACT 


A  finite  element  analysis  was  made  of  crack  grow^th  in  a 
center “ cracked  specimen  subjected  to  monotonically  increasing 
load  until  the  point  of  fast  fracture.  Since  part  of  the 
specimen  experienced  unloading,  the  boundary  value  problem 
which  was  formulated  was  based  upon  incremental  theory  of 
plasticity.  Experimental  load  and  crack-size  records  were 
utilized.  Linear  relations  betw^een  plastic  energy  and  crack 
growth  were  observed.  Fracture  toughness  parameters  G^,  which 
were  evaluated  at  the  onset  of  unstable  crack  propagation,  ob¬ 
tained  from  finite  element  analysis  were  in  good  agreement 
with  those  determined  experimentally. 


1. 


Int  I'oduct  ion 


Linear  elastic  fracture  mechanics  is  widely  used'  by  research 

workers  and  engineers  as  a  tool  to  determine  fracture  toughness, 

K  or  G  ,  for  evaluating  the  susceptibility  of  engineering 
c  c  ’ 

materials  to  unstable  fracture.  Since  ductile  materials  have 
higher  resistance  to  crack  growth  than  brittle  materials,  the 
recent  trend  is  to  develop ’ more  ductile  materials  for  structural 
applications.  However,  ductile  materials  show  a  considerable 
amount  of  crack- tip  plasticity  and  significant  amount  of  crack 
growth  prior  to  the  onset  of  unstable  fracture. 

A  plastic  zone  correction  [1,2],  the  J-integral  [3,4]  and 
COD  methods  [5]  have  been  proposed  to  treat  fracture  involving 
crack-tip  plasticity  in  the  absence  of  crack  growth.  Further, 
the  crack  growth  resistance  curve  method  [6]  has  been  proposed 
for  cases  having  subcritical  crack  growth.  Liebowitz  and 
nftis  [7,8]  introduced  the  nonlinear  energy  method  which  is 
applicable  to  semibrittle  fracture.  Jones  et  al  [9]  developed 
a  corresponding  experimental  program  to  determine  the  toughness 
parameter  and  Liebowitz  et  al  [10]  provided  theoretical 

and  experimental  comparisons  between  their  nonlinear  energy 
method  and  other  existing  methods.  Recently  Jones  et  al  [11] 
used  four  empirical  methods  to  determine  G^  in  the  range  of 
crack  growth. 

In  this  work,  a  finite  element  analysis  is  developed  for 

incremental  theory  of  plasticity. 


crack  growth  problems  based  upon 


In  this  paper,  mention  is  made  of  the  work  of  Newman  [IZ]  whose 
crack  growth  criterion  is  based  on  crack-tip  strain  with  element- 
mesh  size  being  a  parameter,  and  that  of  de  Koning  [15]  who 
found  crack  tip  oj^ening  angle  as  constant  during  the  process  of 
crack  growth.  In  this  analysis,  the  realistic  experimental 
curve  relating  load  and  crack  size  is  taken  as  an  input.  In 
order  to  demonstrate  the  validity  of  this  work,  the  numerical 
output  of  load  and  displacement  is  compared  to  its  experimental 
counterpart.  Moreover,  a  linear  relation  is  observed  between 
plastic  energy  and  crack  size  during  tlie  process  of  crack 
growth.  Tlie  f  inite- element  values  of  ,  based  upon  empirical 
metliods  derived  by  Jones  ct  al  [11],  are  in  good  agreement 
witli  the  experimental  values. 
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2.  Stress-Strain  Relations 


In  linear  elasticity,  the  st ress - strain  relations  for  a 
homogeneous  isotropic  material  can  be  written  in  either  of  the 
following  forms  [14] ; 


“ij  “  ^  bj  ^  • 


(2.1) 


V ^  n  y 


•^ij  "  Til  ^ij  '  ~2u(5X  +  Zy)  kk  ij  ’ 

where  X  and  y  are  the  Lame  constants,  o^.  and  are  the 

stress  tensor  and  strain  tensor  respectively.  Introducing 
stress  deviator  s^j  as  follows: 

""ij  =  ""ij  '■  T  ""kk  '^ij  ’ 


(2.2)  can  be  rewritten  as: 


,  1  -  2v  ^ 

Be.  -  =  (1  +  v)  s-  .  +  —3^ -  0^^  , 


(2.2) 


(2.3) 


(2.4) 


where  B  and  v  are  Young’s  modulus  and  Poisson's  ratio  respec' 


tively . 


In  nonlinear  elasticity,  adopting  the  model  suggested  by 
Ramberg  and  Osgood,  we  generalize  the  stress  -  strain  relations 

as  follows: 


Ee.  .  =  (1  +  V)  s.  -  +  ”  3-—  •Oi.k  '^ij  ^  I  “'"e  ^ij 
where  the  effective  stress  o^  is  defined  as: 

^e  -  I  ^ij  ^ij 


e  ^ij,  (2.5) 


(2.6) 
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In  simple  tension  test,  namely,  all  stress  components  are  zero 
except  =  0,  we  obtain  the  following  nonvanishing  components 

of  c..,  s..  and  the  effective  stress; 

;  .  -‘-J  ■  . 

^  V  ■  ■■  o  ■  :  ^  ' 

he-,  T  =  c  +  ao  , , 


'^^22 


VO  -  2  cto  , 


^11  '5  ^  ’ 


(2.7) 


^22  ^33 


o'e  -  a 


I'or  the  theory  of  incremental  plasticity,  a  certain  portion  of 
the  material  is  said  to  be  in  an  unloading  situation  if  its  current 
effective  stress  0^  is  less  than  the  maximum  effective  stress 
a*  it  experienced  before,  otherwise  it  is  said  to  be  in  a  loading 
situation.  The  relations  between  incremental  stresses  and  in¬ 
cremental  strains  in  the  loading  situation  can  be  derived  from 
eqn.  (2w5)  as  follows; 


=  (1  +  v)  ds..  * 


-  n-1  Q  , 

+  I  aOg  ds. .  +  j  a(n-l)0g  s^-  s^^  ds^^ 


(2.8) 


In  the  unloading  situation,  we  assume  the  incremental  stress- 
strain  relations  arc  the  same  as  those  in  linear  elasticity. 


Bde^j  =  (l+v)ds^^  + 


.do,-,  6.. 

kk  ij 


(2.9) 


The  stress-strain  relation  in  the  case  of  simple  tension  is 
illustrated  graphically  in  Big.  1.  In  the  case  of  generalized 


plane  stress,  equ. 
form: 


-  5  - 

(2.8)  may  be  written  in  the  fallowing  matrix 

I  ’^11  ‘"^12 

"  ^^2  ^^22  ^^23 

'^13  ^23  ^33 

where 

=  {1  +  g  +  h(20^-ay)^}/E  , 

hi2  =  {-  V  -  g/2  +  h(2a,^-ay)  (20^-0^)  }/E  , 

h^3  =  6h(20^-0y)0^y/E  . 

1122  =  ■’■  8  ■*■  ’ 
h^j  =  6h(2ay-o^)o^y/li  , 

=  {2(l+v)  +  3g  +  36h  cJ^y}/E  , 

and 

n-1  n-3 

g  =  ao^  ,  h  =  a(n-l)0g  /4 

Equivalently,  eqn.  (2.10)  can  be  expressed  as 

^11  ^12  "^13 

^12  ^22  ^23 

^13  ^23  ^'33 

where  matrix  [d-^]  is  the  inverse  of  matrix  [h^jj.  In  the  case 
of  unloading,  we  have 
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so  obtained  in  eqns.  (3. 1-3. 2)  are  called  plastic  strains  which 
are  the  residuals  after  the  applied  load  is  removed.  Comparing 
eqns.  ( 3 . 1- 3 . 2)  with  eqn.  ( 2 . 7)  ^  and  eqn.  (2.7)2,  rewrite 

eqn.  (2 . 5)  as 


E.,.  .  (3.3) 

where  the  elastic  strains  e?.  and  plastic  strains  e? .  are 

XJ  J 

Ee®.  =  (l^v)s,.  .  6,.  , 

.5.  n- 1 

Ee?.  =  4-aa  s .  .  . 

ij  2  e  13 

Now  suppose,  in  a  general  case,  the  loading  and  unloading 
history  of  a  certain  portion  of  the  material  can  be  divided  into 
two  stages.  In  the  first  stage  the  effective  stress  increases 
monotonically  to  a*  which  corresponds  to  and  in  the  second 

stage  the  effect  stress  decreases  monotonically  to  which  corre¬ 
sponds  to  Then  tlie  strain  energy  density  cfi  can  be  obtained  as 


o  " 

/ 


-  d  e .  -  + 
1  11 


/ 


a .  -  de . . 
11  11 


_  iri  +  V  ,  1  -  2v,:  ^ 

E^~3—  ""e  ^  —6 - ^ 


n 

— T  oto" 
n  +  1  e 


n+1 


} 


(3.4) 


wc  notice  the  strain  energy  density  can  be  separated  into  two 
parts: 

Irl  +  V  .  1  -  2V,: 

'f’e  =  E^”3~  ""e  ”“6 - ^'^''kk^  ^  ’ 

.n+1 


(3.5) 


L  n  j 

b  =  Try - PT.  ao’- 

p  E(n+1)  e 


(3.6) 
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where  the  elastic  strain  energy  density  4)^  is  a  function  of  the 
current  stresses  S- -  and  the  plastic  strain  energy  density  4)^ 
is  a  function  of  maximum  effective  stress  .  The  irrever¬ 
sible  and  dissipative  nature  of  plastic  strain  energy  density 
is  indicated  in  eqn.  (3.6).  In  two  dimensional  problems,  the  elastic 
and  plastic-strain  energy  per  unit  thickness  over  an  area  A 
can  be  obtained  as: 


/ 

A 


5  /  fe 


/ 


Up  E  I  *p  dA 


(3.7) 


(3.8) 
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4.  Center- Cracked  Specimen 

In  this  work,  we  focus  our  attention  on  a  tectangular 
plate  of  length  2)1,  width  2w,  and  thickness  B,  with  a  centered 
line  crack  of  initial  crack  size  subjected  to  symmetric 

boundary  conditions  (cf.  Fig.  2).  Therefore  only  the  first 
quadrant  of  the  plate  R  =  [x,yl0  <  x  <  w,  0  <  y  5  ^-1  needs  to 
be  analyzed.  Experimentally  one  may  obtain  a  curve  of  applied 
stress  a  vs  crack  size.  A  realistic  experimental  curve  relating 
a  and  a  is  shown  in  Fig.  3.  We  are  interested  in  the  investiga¬ 
tion  of  the  process  that  crack  size  slowly  increases  from  a  -  a^ 
to  a  =  a^  as  the  applied  stress  increases  from  a  =  0  to  a  =  a^. 
Suppose  we  have  two  adjacent  states.  o  -  o  and  a 


State  1,0= 

(2) 
0  ^  ^ 

and 

a  = 

(2) 
a  ^ 

at 

state  2.  Correspondingly, 

the  boundary 

conditions  may  be 

specified  as  (cf.  Fig 

.  4); 

“x 

=  ka 

(i) 

9 

°xy 

=  0 

on 

S^^  =  [X  =  w,  0  <  y 

<  )l]  , 

(4.1) 

“y 

=  ati). 

^xy 

=  0 

on 

S^  s  [y  =  )l,  0  <  X 

<  w]  , 

(4.2) 

u 

=  0 

a 

=  0 

on 

S  =  [X  =  0,  0  <  y 

<  )l]  , 

(4.3) 

X 

9 

xy 

D 

=  0 

9 

°xy 

=  0 

on 

S4  H  [y  =  0,  0  <  X 

<  a^-^M 

,  (4.4) 

=  0 

9 

^xy 

=  0 

on 

S5  E  [y  =  0, 

X  <  w]  . 

(4.5) 

In  other  words,  the  boundary  conditions  associated  with  the 

State  1  to  State  2  can  be  expressed  as  (cf.  Fig.  5) 


process  from 
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6o^  - 

i  «  ^"xy  = 

0 

on 

Si 

y 

(4 

.6) 

0 

(2)-  a 

(1) 

y 

'^xy  = 

0 

on 

Sz 

y 

(4 

.7) 

6u  - 
X 

0 

y 

0 

on 

S3 

y 

(4 

.8) 

60^  - 

0 

,  6a  =  0 

’  xy 

on  r  ^  s 

[y 

= 

0, 

0  <  X  <  a , 

(4 

.9) 

50  = 

0  , 

=  0  on 

[y 

=: 

0,  a(l^<x<a^2) 

], 

(4.10) 

y 

y 

xy 

z 

6u  = 

0 

,  6o 

xy 

0  on  Fj  s 

[y  = 

0, 

a ^2) <x<w]  , 

(4 

.11) 

where  is  the  stress  of  State  1  distributed  along  y  -  0 
between  a^^^<  x  <  a^^\ 

If  we  have  the  solutions  of  State  1,  then,  for  any  given 
adjacent  state  specified  by  ^  and  a  ^  ^ ,  the  boundary  con¬ 
ditions  for  the  transition  from  State  1  to  State  2  are  well  de¬ 
fined  in  cqns.  (4.6-4.11).  The  equations  of  equilibrium  and 
the  strain-displacement  relations  in  the  incremental  form  can 
be  written  as: 


-  0  > 


(4-12) 


(4.13) 


Hqns.  (4.6-4.13)  with  eqn.  (2.8)  and/or  eqn.  (2.9)  mathematically 
define  the  crack  growth  problem  as  a  boundary  value  problem. 

After  the  boundary  value  problem  is  solved,  the  stresses, 
strains,  and  displacements  at  State  2  may  be  obtained  as 


12 


(2)  = 

+  6o •  ■  , 

(4. 

14) 

ij 

ij 

(2)  = 

+  6e  •  •  , 

C4. 

.15) 

ij 

ij 

13 

(2)  = 

+  6u . 

(4 

.16) 

Following  the  same  procedure,  one  may  analyze  stepwise  the 
whole  process  of  slow  crack  growth  up  to  the  onset  of  fast 
crack  propagation. 
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5.  I'inite  l-lcmcnt  Procedure 

Let  a  finite  element  mesh  with  nodal  points  and 
triangular  elements  be  set  up  in  region  R  =  [x,y|0  <  w,  0  <  y  i  ll- 
In  this  work,  N„  =  300.  N  =  518,  and  a  series  of  very  fine 
elements  are  being  distributed  along  the  path  of  slow  crack 
growth.  The  area  pf  those  elements  is  in  the  order  of  10  wS,. 

If  j  is  the  number  of  a  certain  nodal  point,  then  and  U2j 

are  the  displacement  of  that  point  in  x  and  y  direction  respec¬ 
tively,  ^2}  1  ^2j  corresponding  external  concentrated 

force  components  acting  on  that  point.  The  displacement  field 
within  each  triangular  element  is  assumed  to  be  linear  with 
respect  to  the  coordinates.  Therefore  the  strain  field,  and  accord 
ingly  the  stress  field,  within  each  triangular  element  are  constant 
For  each  element  (cf.  Fig.  6),  let  the  incremental  strain  field 
[6el,  incremental  stress  field  [6o]  ,  and  incremental  nodal  point 
displacements  [6]  be  represented  by 


[6el, 

ji 

=  •  '^^y’  “^^xy^  ’ 

(5.1) 

[60] 

T 

=160^,  60y,  6o^y]  , 

(5.2) 

[6]  5 

T 

■  f '^^2i-l’'^^2i  -1  ’'^^2k-l  ’'^“2k^ 

(5.3) 

have 

[6c] 

=  [B][6]  . 

(5.4) 

[6ol 

=  [d] [6e]  , 
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where  [d]  is  a  3  x  3  matrix  as  indicated  in  eqns .  (2.11)  and  (2.12) 

which  are  valid  for  loading  and  unloading  respectively,  and 


[B]  = 


2A 


b.  0  b,  0 


c  . 
1 


c  . 
3 


=i  >^1  b  '’3 


(5.5) 


2A  s  det 


bi  E  Yj  -  y 


1 

1 

1 


k  ’  ^i  ' 


^i  ^i 


y 


'k 


3 

^k 


X .  +  X, 

3  b 


(5.6) 


(5.7) 


with  the  other  coefficients  obtained  by  a  cyclic  permutation 
of  subscripts  in  the  order  of  i,  j,  k.  The  stiffness  matrix 
per  unit  thickness  of  this  triangular  element  linking  incre¬ 
mental  forces  and  incremental  displacements  may  be  obtained  as 

[15,16] : 

T 

[k]  =  [B]  [d][B]  A, 

where  A  is  the  area  of  the  element.  Finally  the  governing 
equation  can  be  written  as; 


2N 


6=1 


K  „  6Uq  =  5f^ 
a6  6  ot 


a 


=  1,2, 


2N 


(5.9) 
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whore  the  2Np  x  stiEfoess  matrix  (K)  is  equal  to  the  sum  of 

N  local  stiffness  matrices  [k]j,  I  -  1 . 

AS  matrix  [J1  of  each  element  depends  on  the  current  stresses 

when  that  element  is  in  a  loading  situation,  the  matrix  [K]  also 
depends  on  the  whole  stress  field.  Thus,  an  iteration  process 
has  to  be  taken  to  solve  the  nonlinear  matrix  equation  (5.9).  The 
iteration  process  could  be  described  as  follows.  Suppose,  in  a  step 
by  step  process,  we  have  the  solutions  of  State  1,  (i.e.,  for  each  ele. 
ment  wo  have  the  current  stresses  and  the  maximum  effective 


tress  o*).  Then,  for  each  element,  we  guess  at  the  stresses  0^^ 
itate  2.  Then  the  current  effective  stress  0^  and  the  effective 

stress  0  maybe  obtained  respectively  from 
c 


%  ■  *  “y  -  "x^y  "  ^“.xy’ 


O.  -  (“v  “v  '  ^°.xy’ 


(5.10) 


(5.11) 


The  corresponding  matrix  (d)  for  each  element  may  be  found  from  one 
of  the  three  following  cases ; 

(1)  =  a*  and  >  a* 

Por  case  (1),  the  element  is  in  a  loading  situation;  matrix  [d] 
assumed  to  be  the  average  matrix  corresponding  to  tiv.o  states  of 

Stresses ,  i.e., 


‘>ij  =  ■ 


(2) 


<  o*  and  f 


(5.12) 
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For  case  (2),  the  element  is  in  an  unloading  situalion,  matrix 
[d]  should  be  the  same  as  indicated  in  eqn.  (2.12). 


(3)  0^  <  and  0^  >  0^ 


For  case  (3) ,  the  clement  has  experienced  unloading  and  is  being 
loaded  again,  the  matrix  [d]  is  obtained  as 


where 


d.  .  =  f  d.  .  (o)  +  (1-f)  ’ 

13  13  ~  ■> 


£  =  (Op  -  Op)/  ■  °e^  ’ 


(5.13) 


(5.14) 


and  d^-(O)  denotes  the  matrix  corresponding  to  vanishing 
stress  state. 

Alter  solving  eqn.  (5.9),  one  may  obtain  the  calculated 
stresses  and  its  corresponding  effective  stress  0^  .  The 

iteration  process  will  be  continued  until,  for  each  elem  , 
the  guessed  stresses  and  calculated  stresses  are  approximate  y 
same  and  the  averaged  percentage  difl.erence  between  0^  and 

is  below  certain  allowable  values  of  error.  After  the  itera¬ 
tion  process  is  completed,  it  is  straightforward  to  calculate 
quantities  of  interest  in  fracture  mechanics. 
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6,  Dissipated  Dnergy  and  Crack  f.rowth 


The  irreversible  and  dissipative  nature  of  plastic  strain 
energy  aensity  have  been  indicated  in  Sec.  3.  After  the  iteration 
process  for  a  certain  incremental  stop  is  completed,  we  have  the 
correctly  gocssed  stresses  Sc.  and  the  corresponding  effective 
stress  of  a  new  state  for  each  finite  element.  If  the  newly 
obtained  effective  stress  is  larger  than  the  maximum  effective 
stress  0*  obtained  in  the  previous  states,  one  should  replace 
O*  by  k  ,  which  will  be  used  as  the  maximum  effective  stress 
for  the  next  incremental  step.  Now  we  recall  the  expiessions 
elastic  and  plastic  strain  energy  density  as  follows: 


♦e  ’  “e  " 


2v 


(o 


kk 


2 

)  )  , 


(6.1) 


•t'p  ' 


n+1 


(6.2) 


Wc  notice  that  (1>^.  depends  on  the  current  stresses  and  therefore 
demonstrates  its  elastic  nature.  On  the  other  hand,  is  a 
monoton  ically  increasing  quantity  because  is  always  increasing, 
this  corresponds  to  the  irreversibility  of  plasticity  and 
the  dissipative  nature  of  plastic  energy.  After  the  finite 
element  analysis  for  each  incremental  step  is  completed,  it 
straightforward  to  obtain  the  total  plastic  energy  P  (per  unit 
thickness)  for  the  whole  region  R  as  follows; 
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N 


P  = 


(6.3) 


i=l 


For  a  stationary  crack  problem,  the  crack  size  is  regarded 
as  a  given  quantity  and  hence  the  question  is  to  determine  the 
critical  load  at  the  onset  of  slow  crack  growth.  However,  as 
the  crack  starts  to  grow  in  a  stable  situation,  one  has  to 
treat  crack  size  as  a  variable.  Strictly  speaking,  one  more 
variable  corresponds  to  one  more  governing  equation  for  the 
system.  Before  one  can  establish  a  governing  equation  for 
crack  size,  the  relation  between  crack  size  and  other  fracture 
parameter  has  to  be  proposed  and  tested  [12,13].  Motivated  by 
the  Griffith’s  concept  [17,18]  that  surface  energy  is  propor¬ 
tional  to  its  area,  and  the  energy  for  the  new  fracture  surface 
has  to  be  supplied  from  the  released  strain  energy  of  the  speci¬ 
men,  we  decided  to  take  the  energy  approach..  Since  crack  growth 
is  generally  regarded  as  a  irreversible  process  and  crack  size 
is  an  ever  increasing  quantity,  we  propose  that  crack  growth 
is  a  function  of  plastic  energy,  i.e.. 


a^  =  f(P) 


(6.4) 


We  mention  that  Broberg  [19]  proposed  a  similar  concept  which 
can  be  expressed  as: 


dUp/da  =  dD^/da  , 


(6.5) 


where  Uq  is  the  energy  flow  to  the  process  region  and  Dq  is 
the  energy  dissipation  in  the  process  region.  However  Broberg  s 
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engineering  approach  is  based  on  the  so  called  approximate 
path- independence  of  J- integral  even  in  the  cases  of  crack 

growth. 

After  wc  analyze  a  realistic  crack  growth  problem  with  all 
the  input  data  obtained  experimentally,  we  observe  the  following 

fact: 

(6.5) 

a  -  =  B(P-l>„) 

which  means  the  amount  of  crack  growth  is  linearl)  proportion 
to  the  increment  of  plastic  energy.  The  normalized  plastic 
energy  per  unit  thickness  p  ^  PE/o^w^  is  plotted  against  normalized 
crack  size  c  5  a/w  in  Fig.  7.  The  correlation  coefficient 
attains  a  high  value  of  0.9997.  Eqn.  (6.5)  may  be  rewritten  as 

p  =  Pq  * 


*  (Pq-  5  “ 


(6.6) 
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7.  Fracture  Toughness 


Liebowitz  and  Eftis  [7,8]  developed  a  nonlinear  energy 
method  to  evaluate  the  toughness  parameter  from  a  single 
load-displacement  record  for  cases  of  no  crack  growth, 
method  may  be  briefly  described  as  follows.  First,  let  the 
experiment  load-displacement  record  be  represented  by 


M 


(7.1) 


where  v  is  the  load  point  displacement,  F  is  the  load,  1/M  is 

the  crack  size  dependent  elastic  compliance,  n  and  k  are  constants. 


Then  for  a  given 


load  F^,  one  may  obtain  a  dimensionless 


quantity  C  representing  the  nonlinearity  as  follows; 


7s  a  .  2nk  f  Csj 

c  -  1  +  ^ 


The  nonlinear  toughness  parameter  is  obtained  as. 


Gc  =  CG^ 


(7.2) 


(7.3) 


where  G  is  the  corre 


spending  linear  elastic  fracture  toughness  and 


may  be  obtained  by  using  Irwin's  G-K  relation,  i 


G  =  K  /E 
c  c 


(7.4) 


And  is 


determined  using  the  following  polynomial  expression 


for  center-cracked  specimen  [20] 


K  =  o  /Fa  [1  -  0*i  ^  (^w^  ^ 

c  ^ 


(7.5) 


Jones  et  al  [11]  further  developed  four  different  methods  to 
account  for  the  effect  of  subcritical  crack  growth.  We  briefly 
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describe  these  four  methods  and  m  each 
experimental  value  of  with  the  numeri 
by  finite-clement  analysis. 


method  we  compare  the 
cal  value  of  obtained 


Method  1 

The  experimental  part  of  this  method  can  be  separated  into 
two  steps:  CD  to  determine  C  from  the  experimental  lead  dis¬ 
placement  record, (2)  to  calculate  by  eqns.  (7. 4-7. 5)  using 

actual  critical  stress  and  initial  crack  size.  For  illus¬ 
trative  purpose,  we  take  specimen  #11  (2024-T3)  as  an  example 
and  the  data  for  that  specimen  are  listed  below; 


2w  =  12  in.  , 

2aQ  =  6  in. , 

1<L  =  32  in.  , 

E  =  10300  ksi, 


(7.6) 


V  =  0.33  , 

=  55.04  ksi  , 
n  =  6.025  , 


n-1 

OtOy 


2  .  207 


And  load-crack  size  relation 


is  plotted  in  Fig-  ^ 


which  is 


taken  as 
data  for 
analysis 


an  input  in  this  work.  As  an  output  the  numerical 
load-displacement  curve  obtained  by  finite  element 
are  plotted  in  Fig.  8  to  compare  with  the  experimental 

In  this  particular  case  we  have 


curve . 


C7.7) 


-  22  - 


(Exp)  =  521  Ib/in.,  (FE)  -  473  Ib/in. 

The  difference  is  about  101;  2%  comes  from  the  calculation  of 
U  due  to  fact  that  eqn.C7.5)  always  gives  higher  value  for 
aL  hence  higher  value  for  from  eqn.  (7.4)  and  8%  comes 
from  the  calculation  of  C  due  to  the  discrepancy  between 
finite-element  and  experimental  data  in  the  load- displacement 

Method  2 


(7.8) 


The  experimental  part  of  this  method  involves  the  deter' 
mination  of  a  hypothetical  load-displacement  curve  assuming 
there  is  no  crack  growth.  Eftis  et  al.  [121  propose 

F  =  M(ag)/M(a^) 

being  used  as  the  critical  load  corresponding  to  the  hypo¬ 
thetical  load-displacement  curve.  Lee  and  Liebowitz  [22]  de 
voloped  a  computer  program  which  can  be  used  to  analyze  a 
center-cracked  specimen  in  the  nonlinear  range  without 
growth,  and  therefore  it  is  straightforward  to  find  F  assuming 
critical  displacement  v^  to  be  the  same  as  obtained  experiment- 

ally.  We  find  that 


G^  (FE)  =  C  (FE)  G^  (FE) 

=  1.1325  X  381.77  (14.01/12.8) 

=  518  Ib/in. 


No  experimental  value  for  G^  has  been  repo 

[11]. 


(7.9) 

rted  by  Jones  et  al 


23 


Method  3 

The  experimental  part  of  this  method  involves  the  deter¬ 
mination  of  My,  ny,  ky,  and  C  accordingly,  from  the  portion  o 
load-displacement  curve,  which  has  no  crack  growth,  and  the 
value  of  F  from  the  following  equation 


f 


^c  M, 


F 


(7.10) 


0 


The  experimental  value  of  is  reported  to  be  520  Ib/in.  and 
the  finite-element  value  of  G  is  found  to  be  519.7  Ib/in. 


Method  4 


Utilizing  k.  and  n  obtained  in  Method  3,  we  recall  eqn. 


.y  -  ..y 

(7.1)  and  write  as  follows 


F,  ’"O 

''c  ■  sf  * 

c  ^ 


(7.11) 


which  is  used  to  determine  M^  M(a^)  .  The  experimental 

part  of  this  method  is  to  determine  C.  based  upon  M^ ,  ny,  ky , 

F^,  and  G^ ,  based  upon  critical  crack  size  a^.  Ihe  experimental 
value  of  G^  is  reported  to  be  581  Ib/in.  However,  we  found 
that  G^  (Fli)  =  491.6  Ib/in.,  which  we  believe  is  more  reliable. 
In  summary,  the  finite-element  values  of  G^  for  four  different 
methods  are  473,  518,  519.7  491.6  comparing  with  experimental 
values  521,  -  ,  520,  581  respectively. 


24 


8.  Conclusions 

In  this  work,  we  focused  our  attention  on  the  finite  ele¬ 
ment  analysis  of  crack  growth  in  a  center-cracked  specimen 
subjected  to  monotonically  increasing  loads  until  fast  fracture 
occurs.  During  the  process  of  crack  extension,  although  the 
applied  load  is  increased  monotonically,  part  of  the  specimen 
has  experienced  unloading.  Therefore,  in  this  work  the  boun  y 
value  problem  formulated  is  based  upon  incremental  theory  of 
plasticity. 

Since,  in  a  crack  growth  problem,  crack  size  is  no  longer 
a  given  constant  parameter  and  must  be  treated  as  a  variable,  a 
governing  equation  for  crack  size  is  needed.  Before  we  can 
formulate  such  a  governing  equation,  an  experimental  curve 
relating  load  and  crack  size  is  taken  as  an  input  and  hopefully 
from  the  output  one  may  find  the  relation  between  crack 
and  other  fracture  parameters  derivable  from  the  stress  field. 
Because  of  the  irreversible  and  dissipative  nature  of  plasticity 
and  crack  growth,  we  believe  that  crack  size  should  be  related 
to  plastic  energy.  Indeed,  after  the  completion  of  the  analysis 
we  do  observe  the  fact  that  the  amount  of  crack  growth  is 
linearly  proportional  to  the  increment  of  plastic  energy.  In 
another  study  of  ours  which  is  in  progress,  the  linear  relation, 
will  be  employed  as  an  input  replacing  the  experimental  load- 
crack  size  curve.  If  the  result  yields  a  load-crack  size 
curve  which  is  in  agreement  with  the  corresponding  experi 


mental,  are  then  one  may  claim  a  governing  equation  for  crack 
size. 

Another  aspect  which  is  of  concern  is  the  fracture  tough¬ 
ness  parameter  .  Liebowitz  and  his  coworkers  [8-10]  developed 
a  nonlinear  energy  method  and  its  corresponding  experimental 
jn'oeedure  for  the  determination  of  at  the  onset  of  stable 
crack  growth.  Recently  Jones  e.t  al  [11]  presented  four  dif¬ 
ferent  empirical  methods  to  account  for  the  effect  of  subcritical 
crack  growth.  These  methods  enable  us  to  determine  the  fracture 
toughness  parameter  G^  at  the  onset  of  unstable  crack  propa¬ 
gation  based  upon  the  experimental  load-displacement  record. 

Then  the  finite- element  values  which  were  determined  for  G^  were 
compared  to  the  corresponding  experimental  values;  satisfactory 
agreement  was  found. 

The  following  points  should  be  stressed.  Inputs  for  the 
computer  program  include  the  material  constants  a  and  n,  in  addi¬ 
tion  to  Young’s  modulus  E  and  Poisson’s  ratio  v,  determined  rrom 
the  simple  tension  test  of  unnotched  specimen.  However,  as  a 
matter  of  fact,  a  and  n  vary  over  a  wide  range  when  one  tries  to 
describe  the  experimental  stress  -  strain  curve  by  a  three  parameter 
expression  and,  from  our  experience,  the  finite  element  computer 
program  is  quite  sensitive  to  the  input  ct  and  n.  We  believe  this 
is  the  major  reason  for  the  existing  discrepancy  between  the  exper 
mental  load-displacement  curve  and  its  finite-element  counter¬ 
part,  which  has  a  5.2%  difference  in  displacement  at  the  critical 
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load.  Although  tlicrc  is  a  discrepancy,  we  find  that  the  finite- 
element  values  of  based  on  the  empirical  methods  have  less 
scatter  than  the  experimental  counterparts.  Generally  speaking, 
this  work  provides  support  to  the  nonlinear  energy  method  de¬ 
veloped  by  Liebowitz  and  hftis  [7],  Eltis  and  Liebowitz  [8]  and 
its  corresponding  experimental  procedures  Jones  et  al  [9]  , 
Liebowitz  et  al  [10]  and  Jones  ct  al  [11],  in  determining  the 
fracture  toughness  values  in  the  range  of  stable  crack  growth. 
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Fig.  3  Load  (Stress)  vs  Crack  Size;  2W=12  in.,  2L='32  in,, 
B=0.062  in.,  E-10300  ksi,  0^=55.04  ksi,  n=6.025, 

v=0.33,  aoY'’'‘'=2.207;  k=0 


Fig.  6  Coordinates,  Incremental  Displacements 
and  Forces  of  the  j-th  Nodal  Point 
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Material  :  2024-T3 
Correlation  Coefficient  =0.9997 


Fig.  7  Plastic  Strain  Energy  vs  Crack  Size:  2w-12  in.,  2^32  in.,  2ao  6  in., 

8=0.062  in.,  E=10300  ksi,  v=0.33,  n=6.025,  Oy=55.04  ksi,  ttOy  "  =2.207:  k=0 


0  2  4  6  8 


Displacement  (10"*  in.) 

Fig.  8  Load-Displacement  Curve,  2w=12  in., 
2L=32  in.,  2ao=6  in.,  B=0.062  in., 
£=10300  ksi,  Oy=55.04  ksi,  v=0.33, 

n=6.025.aoy""^  =2.207;  k=0 
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